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Abstract. Precision studies of QCD at e + e~ colliders are based on measurements of event 
shapes and jet rates. To match the high experimental accuracy, theoretical predictions to next- 
to-next-to-leading order (NNLO) in QCD are needed for a reliable interpretation of the data, 
l' We report the first calculation of NNLO corrections {0(a 3 s )) to three-jet production and related 

D ' event shapes, and discuss their phenomenological impact. 



1. Introduction 

Measurements at LEP and at earlier e + e~ colliders have helped to establish QCD as the theory 
of strong interactions by directly observing gluon radiation through three-jet production events. 
The LEP measurements of three-jet production and related event shape observables are of a very 
high statistical precision. The extraction of a s from these data sets relies on a comparison of the 
t^- \ data with theoretical predictions. Comparing the different sources of error in this extraction, one 

finds that the purely experimental error is negligible compared to the theoretical uncertainty. 
There are two sources of theoretical uncertainty: the theoretical description of the parton-to- 
hadron transition (hadronisation uncertainty) and the uncertainty stemming from the truncation 
of the perturbative series at a certain order, as estimated by scale variations (perturbative or 
scale uncertainty). Although the precise size of the hadronisation uncertainty is debatable 
and perhaps often underestimated, it is certainly appropriate to consider the scale uncertainty 
as the dominant source of theoretical error on the precise determination of a s from three-jet 
observables. 

The three-jet rate and event shapes related to it can be expressed in perturbative QCD by 
dimensionless coefficients. These coefficients depend either on the jet resolution parameter or on 
the event shape variable. Typically, one denotes these coefficients by A, B,C, . . . at LO, NLO, 
NNLO, etc. 

The perturbative expansion for the distribution of a generic observable O up to NNLO for 
renormalisation scale /i 2 = s and a s = a s (s) is given by 
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Figure 1. Results of a NLO QCD fit to 18 hadronic observables for (a) a fixed renormalisation 
scale \x = ->/i and (b) optimising the scale choice to fit the data. The figures are taken from 
Ref. [1] 



normalised to the tree-level cross section for e + e — > qq. It can easily be related to the physical 
observable 4^ using the known relation between <7h a d an d <?o- 

Ignoring the numerically negligible singlet contribution, A, B and C depend only on the jet 
resolution parameter or the event shape variable under consideration, and are independent of 
electroweak couplings, centre-of-mass energy and renormalisation scale. 

QCD studies of event shape observables at LEP pQ based around the use of fixed-order NLO 
parton-level event generator programs [2] have shown that the current error on a s from these 
observables [3] is dominated by the theoretical uncertainty. To illustrate this, the DELPHI 
collaboration made fits to 18 hadronic observables (a) at fixed renormalisation scale, [i = yfs 
and (b) by treating the renormalisation scale as a free parameter in the fit. The results are 
displayed in Fig. [H taken from Ref. [1] . With a fixed scale (fig [U(a) ) , the value of the strong 
coupling extracted varies considerably amongst leading to a value of a s (Mz) = 0.1232 ± 0.0116. 
On the other hand, fig H!(b) shows the result of letting the scale vary (and thereby estimating 
the uncalculated higher order corrections), and leads to a more consistent fit and a much smaller 
error a s {Mz) = 0.1168 ±0.0026. Clearly, to improve the determination of a s , the calculation of 
the NNLO corrections to these observables becomes mandatory. 

2. Thrust distribution 

As an example, we focus on the thrust distribution, O = T. The LO and NLO coefficients 
A(T) and B(T) are displayed for comparison in Figure [2j In the numerical evaluation, we use 
M z = 91.1876 GeV and a s (M z ) = 0.1189 [3]. 

The NNLO coefficient C (T) has been recently computed |5] . It is based on two-loop 7* — > qqg 
matrix elements [6j [7] , one-loop four-parton matrix elements [8] and tree- level five-parton matrix 
elements [9]. 

The two-loop 7* — * qqg matrix elements were derived in [6] by reducing all relevant 
Feynman integrals to a small set of master integrals using integration-by-parts [10] and Lorentz 
invariance [TT] identities. The master integrals [12] were computed from their differential 
equations [H] and expressed analytically in terms of one- and two-dimensional harmonic 




Figure 2. Coefficients of the leading order, next-to- leading order and next-to-next-to- leading 
order contributions to the thrust distribution. 



polylogarithms [T3] . 

The one-loop four-parton matrix elements relevant here [8] were originally derived in the 
context of NLO corrections to four-jet production and related event shapes [HI [15] . 

The four-parton and five-parton contributions to three-jet-like final states at NNLO contain 
infrared real radiation singularities, which have to be extracted and combined with the infrared 
singularities present in the virtual three-parton and four-parton contributions to yield a finite 
result. In our case, this is accomplished by introducing antenna subtraction functions [T6| [T7], 
which encapsulate all singular limits due to the emission of one or two unresolved partons 
between two colour-connected hard partons, and are sufficiently simple to be integrated 
analytically [18j . 

The resulting numerical programme, EERAD3, yields the full kinematical information on a 
given multi-parton final state. It can thus be used to compute any infrared-safe observable 
related to three-particle final states at O(o%). The NNLO coefficient C(T) for the thrust 
distribution obtained with EERAD3 is shown in Fig. [5j. 

Figure O displays the perturbative expression for the thrust distribution at LO, NLO and 
NNLO, evaluated for LEP and ILC energies. The error band indicates the variation of the 
prediction under shifts of the renormalisation scale in the range \i 6 [Q/2;2Q] around the e + e~ 
centre-of-mass energy Q. 

It can be seen that even at linear collider energies, inclusion of the NNLO corrections enhances 
the thrust distribution by around 10% over the range 0.03 < (1 — T) < 0.33, where relative scale 
uncertainty is reduced by about 30% between NLO and NNLO. Outside this range, one does 
not expect the perturbative fixed-order prediction to yield reliable results. For (1 — T) — > 0, 
the convergence of the perturbative series is spoilt by powers of logarithms ln(l — T) appearing 
in higher perturbative orders, thus necessitating an all-order resummation of these logarithmic 
terms [191 [20] , and a matching of fixed-order and resummed predictions [21] . 

The perturbative parton-level prediction is compared with the hadron-level data from the 
ALEPH collaboration [22] in Figure The shape and normalisation of the parton level NNLO 
prediction agrees better with the data than at NLO. We also see that the NNLO corrections 
account for approximately half of the difference between the parton level NLO prediction and 
the data. 

3. Conclusions 

We developed a numerical programme which can compute any infrared-safe observable through 
to 0(o!g), which we applied here to determine the NNLO corrections to the thrust distribution. 
These corrections are moderate, indicating the convergence of the perturbative expansion. Their 
inclusion results in a considerable reduction of the theoretical error on the thrust distribution 
and will allow a significantly improved determination of the strong coupling constant from jet 
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Figure 3. Thrust distribution at LEP and at the ILC with Q = 500 GeV. 



observables from existing LEP data. 
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